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With future galaxy surveys, a huge number of Fourier modes of the distribution of the large 
scale structures in the Universe will become available. These modes are complementary to those 
of the CMB and can be used to set constraints on models of the early universe, such as inflation. 
Using a MCMC analysis, we compare the power of the CMB with that of the combination of 
CMB and galaxy survey data, to constrain the power spectrum of primordial fluctuations generated 
during inflation. We base our analysis on the Planck satellite and a spectroscopic redshift survey 
with configuration parameters close to those of the Euclid mission as examples. We first consider 
models of slow-roll inflation, and show that the inclusion of large scale structure data improves the 
constraints by nearly halving the error bars on the scalar spectral index and its running. If we 
attempt to reconstruct the inflationary single-field potential, a similar conclusion can be reached 
on the parameters characterizing the potential. We then study models with features in the power 
spectrum. In particular, we consider ringing features produced by a break in the potential and 
oscillations such as in axion monodromy. Adding large scale structures improves the constraints on 
features by more than a factor of two. In axion monodromy we show that there are oscillations with 
small amplitude and frequency in momentum space that are undetected by CMB alone but can be 
measured by including galaxy surveys in the analysis. 



I. INTRODUCTION 

Inflation [U-HI is a successful paradigm to explain 
the observed Cosmic Microwave Background (CMB) 
anisotropies. Currently, the Planck satellite Q is taking 
data and one of its main goals will be to accurately mea- 
sure the primordial power spectrum of scalar and tensor 
fluctuations, with the aim of constraining the inflationary 
parameters and, possibly, ruling out models. Indeed, the 
spectrum of primordial perturbations represents an im- 
portant source of information on the inflationary phase: 
during inflation scalar and tensor fluctuations are pro- 
duced with amplitudes and shapes related to the dynam- 
ics of the fields driving inflation. For instance, slow-roll 
inflation predicts a firm deviation from a purely scale- 
invariant spectrum, i.e. n s = 1. Indeed, by combining the 
WMAP CMB data with BAO distance measure @ and 
the Hubble constant Ho measurement @, the WMAP 
team reported a deviation of n s — 1 = — 0.037 ± 0.012, a 
measurement excluding the purely scale-invariant spec- 
trum by more than 3cr Q. 

On the other hand, a series of new surveys is being 
planned to accurately measure the large scale structure 
of the Universe. The number of Fourier modes avail- 
able in these surveys will be larger than those collected 
by the Planck satellite. Furthermore, these modes will 
be on scales smaller than those probed by the CMB (on 
scales where primary anisotropies dominate) and sample 
the full three dimensional (3-D) structure of the den- 
sity field. Even with the drawback that the late-time 
density field has evolved gravitationally and that non- 
linearities severely limit the amount of information that 
can be extracted from large-scale structure surveys, their 



large number of modes make them an invaluable observ- 
able. Hence, galaxy surveys are complementary probes 
to the CMB anisotropies and any attempt of constraining 
the initial conditions should take them into account. The 
literature on adding large-scale structure to CMB data 
to constrain inflationary model is extensive, e.g., [§14121] 
and references therein. 

In this paper we perform a Markov Chain Monte Carlo 
(MCMC) analysis of the constraining power on the pri- 
mordial power spectrum gained by combining Planck 
data with those from a future spectroscopic galaxy red- 
shift survey with Euclid-like [13| characteristics. This 
approach is more time-consuming and numerically inten- 
sive than the popular Fisher matrix approach to fore- 
casts, but it is more precise and robust and solves many 
of the drawbacks of Fisher-based forecasts. As explained 
in Sec. [H] for Planck we use 2.5 years multiple chan- 
nel mock data fl4| while for the galaxy survey we use 
a survey with specifications described in Sec III Bl as an 
example of stage IV surveys, according to the classifica- 
tion of [IH . Such specifications are similar (although not 
identical) to those of Euclid [H, [TtJ ■ Hereafter we refer 
to such a data set as LSS. 

We only consider single-field models of inflation where 
the inflaton field which drives inflation is a canonical 
scalar field characterized by some potential. In the liter- 
ature there are mainly three approaches to constraining 
inflation. The simplest is to assume that the inflaton 
field underwent a period of slow-roll during the phase 
where the observable primordial perturbations were gen- 
erated, as reviewed in Sec. IIII Al This allows to connect 
observations-based parameters, which directly character- 
ize the primordial power spectrum, such as its amplitude, 
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its tilt and the running of the tilt, to potential-based pa- 
rameters, written in terms of slow-roll quantities. We 
follow this approach in Sec. IIIIB| where we show that 
LSS surveys improve the constraints on both the spec- 
tral tilt and its running by a factor of about 1.5. 

A more sophisticated approach is to try to reconstruct 
the inflaton potential from the data. Many methods, re- 
viewed in Sec. lIII Cl have been proposed in the literature. 
All of them require some assumptions on the smoothness 
of the potential and their conclusions always depend on 
these priors. In our potential reconstruction we will as- 
sume that the inflaton potential remains smooth over a 
field range of the order of the Planck mass and, following 
Refs. [laQjl, that also the first two slow- roll parameters 
- those already constrained by the data - remain small. 
We use this method to forecast the power of reconstruct- 
ing the inflaton potential up to its fourth derivative. For 
the first three derivatives we reach analogous conclusions 
to those drawn in the slow-roll case, although worsened 
by the addition of one extra parameter. 

The third approach is to assume a specific model - 
given in terms of the inflaton potential - characterized by 
a handful of parameters. In this approach one estimates 
the values and errors of these parameters from the data. 
We use this scheme in Sec.|lV]to study models generating 
features in the primordial spectrum. Local signatures in 
the power spectrum can arise, for instance, when the in- 
flaton potential has sharp features [20T - I22T ] , when there is 
a transition between different stages in the inflaton evolu- 
tion [23|, Hi| , when more than one field is present 2f| H(| , 
from particle production during inflation [13, [H | , mod- 
ulated preheating [2!|[3(|, or, more recently, in models 
motivated by monodromy in the extra dimensions (3l| 
(see also |32()- These features represent an important 
window on new physics because they are often related 
to UV scale phenomena inaccessible to Earth-based ex- 
periments. Furthermore, as single-field slow-roll inflation 
predicts a smooth power spectrum, it is important to es- 
timate our power of falsifying this simplest model [33l - l37| . 

In Sec. II V Al we discuss what are the main theoretical 
limitations to measure features using the CMB and the 
large-scale structure. For the large-scale structure these 
are the size of the volume of the survey, which limits the 
smallest measurable width of the features in momentum 
space, and the number of modes contained in the survey. 
The latter is limited by our understanding and control 
of the small non-linear scales, where a large number of 
modes is present. The CMB angular power spectrum is 
an integrated observable: the Ci are a 2-D quantity aris- 
ing from a convolution of the 3-D primordial power spec- 
trum with the radiation transfer function. Small scales 
features get smoothed by this projection process and may 
remain unresolved [H}. 

Two typical examples of features in the primordial 
power spectrum are ringing signatures from a step in 
the potential [2(| (Sec. IIVB[) and superimposed oscilla- 
tions from axion monodromy HH [H, E3| fSec. UVCl) . Wc 
discuss these two cases and find that adding large-scale 



structure to the CMB may be crucial in detecting fea- 
tures. In particular, in the case of oscillations we show 
that, as a function of the frequency, the likelihood is very 
different from a Gaussian. Thus, detecting the correct 
oscillation becomes more similar to the process of tun- 
ing a radio channel than to a Fisher matrix analysis. In 
this case we show that there are features remaining un- 
detected even by variance-limited CMB data that can be 
measured with very high precision when LSS data are 
added. In Sec. [V] we draw our conclusions. 

For the MCMC analysis employed in this paper we 
used a modified version of CosmoMC [ijj . Contrarily to 
the original code, this modified version allows some of the 
parameters to have periodic boundary conditions. More- 
over, for an accurate computation of the CMB angular 
power spectrum from a primordial spectrum with fea- 
tures, the standard publicly available CMB codes j42Tj45j 
are not sufficient. One of us (ZH) has developed a new 
numerical CMB integrator, built ad hoc to treat sharp 
features and frequent oscillations. The whole package is 
described shortly in Appendix [X] and more extensively in 
Ref. 0. 



II. FORECASTS FORMALISM AND 
OBSERVABLES 



In the following we perform a MCMC calculation of 
the 68.3% and 95.4% confidence level (CL) forecasts, us- 
ing "mock" power spectra and the associated predicted 
errors for the surveys considered. For the CMB we use 
the expected signal to noise of the Planck satellite [l4| 
while, as an example of upcoming galaxy surveys, we 
focus on a stage IV galaxy survey (following the conven- 
tion of Ref. [la]). In particular, as "straw- man" survey 
we use the specifications reported in Sec. Ill Bl which are 
similar to those of Euclid jl6l . Ht| . However, our find- 
ings will qualitatively apply to this more general class 
of surveys, which are likely to have comparable perfor- 
mances. Furthermore, we only consider the spectroscopic 
component of the survey and only make use of the galaxy 
power spectrum. In particular, we will not directly use 
the weak gravitational lensing power spectrum, although 
weak lensing measurements will be useful in confirming 
some of the assumptions used here, such as for instance 
the scale-independence of galaxy bias on linear scales. 

As fiducial cosmology for our analysis we use a flat 
ACDM model with Q c h 2 = 0.1128, Q b h 2 = 0.022, 
h = 0.72, cr 8 = 0.8 and r rc = 0.09, where fit, and f2 c 
are the energy fraction of baryons and cold dark matter, 
respectively, at redshift zero, h is the Hubble constant in 
units of 100 km/(s Mpc) and r re is the reionization optical 
depth. 
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TABLE I. Planck Instrument Characteristics for 30 months 
of integration. 



Channel Frequency (GHz) 


70 


100 


143 


Resolution (arcmin) 


14 


10 


7.1 


Sensitivity - intensity (^K) 


8.8 


4.7 


4.1 


Sensitivity - polarization (fiK) 


12.1 


7.5 


7.8 



A. Forecasts for Planck 

For our forecast analysis we use the expected signal to 
noise of Planck 2.5 years (5-sky surveys) of CMB multiple 
channel data. We use three channels for Planck mock 
data and we assume that the other channels are used for 
foreground removal and thus do not provide cosmological 
information. We list the instrument characteristics of the 
channels used in our analysis in Table HI We take the 
detector sensitivities and the values of the full width half 
maxima from the Planck "Blue Book" [ljj] ■ 

Given a likelihood function £, we define \ 2 = — 21n£. 
For a nearly full-sky CMB experiment (here we use an 
observed fraction of sky / s k y = 0.75) x 2 can be approxi- 
mated by 0, 0] 



]T (2e + i)f sky 



nBB 



C 



BB 



In 



C 



BB 
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'cJ T c EE ~{Cj E ) 2 \ 
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nTTnEE 
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where we assume / m j n = 3 (given the sky cut the error 
on I = 2 is large) and ^ max = 2500. In this formula, 
Cf Y are the model-dependent theoretical angular power 
spectra for the temperature, E and B polarizations and 
their cross-correlations, with AT, Y = {T, E, B}. They are 
given by Cf Y = Cf Y +NjF Y , where the d are calculated 
using the publicly available code CAMB [42| and the N# 
are noise spectra, which we compute assuming Gaussian 
beams. The estimators of the measured angular power 
spectra, Cf Y , include the contribution from the noise. 

We use the model introduced in [47[ (and later updated 
in [1§| ) to propagate the effect of polarization foreground 
residuals into the estimated uncertainties on the cosmo- 
logical parameters. For simplicity, in our simulation we 
consider only the dominant components in the frequency 
bands that we arc using, i.e., the synchrotron and dust 
signals. We assume that foreground subtraction can be 
done correctly down to a level of 5%. 



B. Forecasts for LSS surveys 

We now turn to the LSS survey forecasts. 



We model the galaxy power spectrum in redshift space 
as (e.g., HUll) 

P g {k, n; z)=(b+ J> 2 ) 2 D 2 {z)P m {k) cxp (-k 2 fj, 2 a 2 ) , 

(2) 

where /i is the cosine of the angle between the wavenum- 
ber k and the line of sight, D(z) is the linear growth 
factor, / = d\nD/dlna is the linear growth rate, P m (k) 
is the matter power spectrum today (at redshift z = 0) 
and <j r parameterizes the effect of small scales velocity 
dispersion and redshift errors as explained below. 

The term f/j, 2 accounts for the redshift distortions due 
to the large-scale peculiar velocity field [49) . which is 
correlated with the matter density field. The exponen- 
tial factor on the right-hand side accounts for the ra- 
dial smearing due to the redshift distortions that are un- 
corrected with the large-scale structure. In particular, 
we consider two contributions. The first is due to the 
redshift uncertainty of the spectroscopic galaxy samples 
which is estimated to be a z = 0.001(1 + z) [l6| . The sec- 
ond comes from the Doppler shift due to the virialized 
motion of galaxies within clusters, which typically has a 
pairwise velocity dispersion a g of the order of few hun- 
dred km/s. This can be parameterized as ^(1 + z) [5l| . 
Taking the qudratic mean of the two contributions, this 
turns into a comoving distance dispersion a r given by 



(1 + ^ 
H 2 (z) 



(10- 6 + a 2 /2) 



(3) 



where H is the Hubble parameter as a function of the red- 
shift. Note that the two quadratic velocity dispersions in 
the parenthesis are of the same magnitude. Practically, 
neither the redshift measurement nor the virialized mo- 
tion of galaxies can be precisely modeled. In particular, 
the radial smearing due to peculiar velocity is not neces- 
sarily close to Gaussian. Thus, eq. ((2]) should not be used 
for wavenumbers k > H ,\ z ) \ , where the radial smearing 
effect is important. This consideration strictly applies 
only to the radial direction. However, gravitational non- 
linearities affect equally the radial as well as the angular 
direction. By imposing this cut at all fi we make sure 
that we exclude from our analysis non-linear scales. We 
thus introduce a UV cutoff fc max as the smallest value be- 
tween — — r and 7^5, where R is chosen such that the 
r.m.s. linear density fluctuation of the matter field in a 
sphere with radius R is 0.5. The resulting maximum wave 
numbers /c max arc reported in Table [H] for each redshift 
bin defined below. Note that these UV cutoffs are rela- 
tively conservative, even comparing to the current values 
used for the SDSS catalog at redshift ~ 0.35 [52| . 

The survey volume is split into 8 redshift bins. The 
redshift ranges and expected numbers of observed galax- 
ies per unit volume are given in Table|ITJ For the observed 
sky we use 20'000 sq. degrees (Euclid is expected to only 
cover 15'000 sq. degrees, so this represents a departure 
from the baseline for Euclid) . With current specifications 
the Euclid mission will not include the first redshift bin; 
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TABLE II. Redshift bins used in the analysis. With current specifications, the Euclid mission will not include the first redshift 
bin; other ground-based surveys will, albeit with smaller sky coverage. 



z 


Az 


n obs (/i 3 Mpc" 3 ) 


fc min (/iMpc" 1 ) 


fe ma x (/iMpc" 1 ) 


bias b 


0.6 


0.2 


3.56 x 10~ 3 


0.0033 


0.15 


1.053 


0.8 


0.2 


2.42 x 10~ 3 


0.0029 


0.17 


1.125 


1.0 


0.2 


1.81 x 10~ 3 


0.0027 


0.20 


1.126 


1.2 


0.2 


1.44 x 10~ 3 


0.0026 


0.21 


1.243 


1.4 


0.2 


0.99 x 10~ 3 


0.0025 


0.22 


1.292 


1.6 


0.2 


0.55 x 10~ 3 


0.0024 


0.22 


1.497 


1.8 


0.2 


0.29 x 10~ 3 


0.0024 


0.23 


1.491 


2.0 


0.2 


0.15 x 10~ 3 


0.0024 


0.23 


1.568 



other -ground based- surveys will, albeit with smaller 
sky coverage: ~ lO'OOO sq. degrees. However, the dif- 
ference in total survey volume between assuming lO'OOO 
sq. degrees or 20/000 sq. degrees in the first redshift bin is 
only few percent and this has virtually no impact on our 
results. The number density of galaxies that can be used 
is n = e ■ n bs, where e is the fraction of galaxies with 
measured redshift. Due to the high accuracy of the spec- 
troscopic redshift and the width of the bins, we ignore 
the bin-to-bin correlations and write x 2 as 



x 2 = 



E 

k,fi,z bins 



Pg, model Q-% fiducial 
^Pg, fiducial , 



(4) 



As on large scales the matter density field has, to a 
very good approximation, Gaussian statistics and un- 
corrected Fourier modes, the band-power uncertainty is 
given by [El [13 



2(2tt) s 



(2TTk 2 dkdfi)(Anr 2 f sk ydr) 



1/2 



1 

P 9 +- 

n 



(5) 



where r is the comoving distance given, for a flat FRW 
universe, by r(z) = J Q cdz' / H(z'). The second term in 
the parenthesis is due to shot noise, under the assumption 
that the positions of the observed galaxies are generated 
by a random Poisson point process. In practice n is not 
known a priori and is calibrated by galaxies themselves. 
The imperfect knowledge of n can bias P g on the scale 
of the survey [53j]. This is taken into account by using 
an IR cutoff fc m i n ~ Gpc -1 , chosen such that fc m i n (z) = 
27r/V 1 / 3 (z), where V(z) is the comoving volume of the 
redshift slice z — Az/2 < z < z + Az/2. If not otherwise 
specified, in each redshift bin we use 30 fc-bins uniformly 
in In k and 20 uniform /i-bins. For models with features in 
the primordial power spectrum, we use a different binning 
scheme in fe, as explained in Section ITVl For each redshift 
bin the value of fe m i n is reported in Table HT1 

We consider a pessimistic and an optimistic case. For 
the pessimistic case we take e = 0.35. For the fiducial 
value of the bias, in each redshift bin we take the one 
reported in the last column of Table [TTJ We assume that 
o~g is redshift dependent and choose <r g = 400km/s as 



the fiducial value in each redshift bin 55|. Then, we 
marginalize over b and a g in the 8 redshift bins, for a 
total of 16 nuisance parameters. For the optimistic case 
we take e = 0.5. For the fiducial value of the bias we 
use b = \/l + z, which is about 10-30% larger than the 
one used in the pessimistic case. In this case we assume 
that <j g is redshift independent, so that there are only 9 
nuisance parameters, i.e. b in 8 redshift bins and o~ g . 

Gaussian statistics for the matter density field and a 
scale-indepedent bias in each redshift bin should be good 
approximations. On the other hand, we have discarded 
all information beyond the linear scales and marginalized 
over the bias. In reality, weak lensing could allow us to 
get a good prior knowledge on the galaxy bias. Moreover, 
it is likely that we can obtain more information by mod- 
eling the mildly non-linear regime, by taking into account 
the bin-to-bin cross correlations (tomography) and using 
the three-point correlations. These topics are beyond the 
scope of this paper. In conclusion, the reader should bear 
in mind that the forecast made here is "clean" (linear- 
scale only) and hence "conservative" [56| . 



III. CONSTRAINING SMOOTH 
INFLATIONARY POTENTIALS 

A. Primordial perturbations from inflation 

Let us start by introducing the standard parameteri- 
zation of the primordial perturbations after inflation (see 
for instance j57| )■ In a perturbed flat FRW universe, by 
choosing a gauge where there are no fluctuations in the 
inflaton field we can write the spatial part of the metric 
as 

gij = a 2 (t) [(1 + 2C) Sij + hij] , ha = , <9;/iy = , 

(6) 

where a is the scale factor, ((t, x) is a scalar perturbation 
which has the property to be conserved on super-horizon 
scales for adiabatic perturbations, and hij (t, x) is a spin- 
2 quantity charactering the tensor modes of the metric 
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The power spectrum of primordial scalar perturba- 
tions, P((k), is defined by 



(CkCk<> =(2TT) 3 S D (k + k / )P C (k) , 



(7) 



where (...) denotes an ensemble average. It is useful to 
define a dimensionless spectrum for scalar fluctuations, 



V s {k) = ^P<(k) 



(8) 



The deviation from scale-invariance of the scalar spec- 
trum is characterized by the spectral index n s , defined 
by 



f 



dlnV s 
G?ln k 



(9) 



where n s = 1 denotes a purely scale-invariant spectrum. 
We also define the running of the spectral index a s as 



dn s 
din k 



(10) 



These quantities are defined at a particular pivot scale, 
which for our analysis we chose to be fc* = 0.05 Mpc -1 . 
Thus, with these definitions the dimensionless power 
spectrum can be written as 



V s {k) = A s 



s (fc„)-l+ia 3 (fc«)ln(fc//c„) 



(11) 



where A s is a normalization parameterizing the ampli- 
tude of the fluctuations. 

The tensor quantity hij contains two modes of different 
polarization, h + and h x , each with power spectrum given 
by 



(h k h k ,) = (2Tr) 3 5(k + k')P h (k) 



(12) 



Defining the dimensionless power spectrum of tensor fluc- 
tuations as Vtik) = 2^2-P/ l (fc) (the factor of 2 comes 
from the two polarizations), the ratio of tensor-to-scalar 
fluctuations is given by 



(13) 



During inflation the fractional change of the Hub- 
ble rate per e-fold is small, e = —H/H 2 <C 1, where 
H = dhia/dt is the Hubble rate and a dot denotes 
derivative with respect to the cosmic time. Thus, the 
above observables take a simple form at leading order in 
e. The spectrum of fluctuations is given by (see [5^] for 
an expression beyond leading order in e) 



V s (k) 



1 H 2 



8n 2 e M 2 



(14) 



,H 



where M p = (87rG) -1 / 2 is the reduced Planck mass 
and the right-hand side is evaluated when the comov- 
ing scale k exits the Hubble radius, while the spectral 



index ((9]) is given, using eqs. (|9]) and (|T4|) . by n s = 
1 — 2e — t/(He)\k =a H- The relative variation of e in a 
Hubble time is typically small, so that the power spec- 
trum is close to scale invariance. The spectrum of tensor 
fluctuations is given by 



9 H 2 



(15) 



iH 



which shows that the ratio of tensor-to-scalar fluctua- 
tions in eq. (fT3")) is simply related to the first slow-roll 
parameter by r = 16e. 

In slow-roll inflation, during the adiabatic evolution 
one can related e and e to the potential-based parameters 
ey and rjv-i defined as 



M 2 



9 V" 
M — 

p V 



(16) 



where a prime denote the derivative with respect to the 
inflaton field. In terms of these parameters, the scalar 
spectral index reads, at leading order in slow-roll, 



n s = 1 — 6ey + 2rj\ 



(17) 



Another observational parameter is the running of the 
spectral index introduced in eq. ([TU|. which at leading 
order in slow-roll can be rewritten as 

a s = 16ej] V - 2Ae 2 v - 2£ v , (18) 

where £y is the third slow-roll parameter, defined as 

„ V'V" 

tv = Af p 4 -^ . (19) 
The running is thus second-order in slow-roll. 



B. Scalar spectral index and running in slow-roll 
models 



Let us compute the forecast on the spectral index n s 
and the scalar running a s from slow-roll inflation. As 
a fiducial model we consider chaotic inflation based on 
the quadratic inflaton potential V = ^m 2 (f> 2 [6(J. In this 
case eq. (JIBJ gives e v = tj v = 2M 2 /4> 2 = 1/(2 A) (while 
from eq. (fT9")l £y = 0) and at leading order in slow-roll 
one finds, using eqs. (fT5"|) . (TT7)) and (TT5)) . 

n s = l-2/A, r = 8/N, a s = -2/N 2 , (20) 

where A is number of e-folds from Hubble crossing to 
the end of inflation. Choosing A = 62.5 this yields n s = 
0.968, r = 0.128 and a s = as our fiducial model. This 
choice is fully consistent with a joint analysis of the latest 
cosmological data @, [H, HH, Ull • 

We consider the forecast constraints on eight cosmo- 
logical parameters, i.e. fifc/i 2 , Q c h 2 , 8, r rc , ln^4 s , n s , a s , 
and r. Here 9 is the angle extended by sound horizon 
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TABLE III. Cosmological Parameters 





Planck 


Planck + LSS opt. 


Planck + LSS pess. 


Planck + LSS pess. Fisher 


9, b h 2 


u - uzzul -0. 00012 


, , io nnoos 
U.UZZUU_ 00008 


, , io nnoos 
f)2900 + 
U.UZZUU_ 00008 


, io OOOOS 
f)29n0 + 
u.uzzuu_ 00008 


Q c h 2 


0.1127±g;ggSg 


0.11280lg;™023 


_ . . _ „ . I n 0009H 

U.llZOl_ 00026 


„ . . ___i_n nniWi 
U.11ZO_ .00025 





., ^^^0 0009 

1 04fi0 + 
i.uy:OU_ 0002 


., „,*-l_0 0009 

1 04fi0 + 

l.U4:OU_ QQQ2 


., „,*-l_0 0009 

1 04fi0 + 

l.U4:OU_ QQQ2 


., _..in OOOO 

1 04fi0 + 
i.uy:OU_ 0002 


Txe 


0SQQ + 

U.UOTO_q .0038 


u.uyuu_ 0030 


„ „„„„_i_n oo'^o 
U.UJUU_ 0030 


u.uyuu_ 0023 


n s 


0.96821^029 


_ _ _ _ . j_n on 1 7 

0.9681±g;ggi? 


0-968lHo2o 


_ _ „„_|_0 001 s 

0-968lHoi8 




-0.00051^056 


-0.000±°;°030 


„ _ _i_o 00^9 


^ /-innn-t-0 0.02R 

o.ooooin° 2 » 


ln(10 10 A) 


o nlq +o.oo9 

°- ul3 -0 .008 


o Olq +0.006 


o.ury_o.oo6 


o nl q+o oo5 
o.uiy_ 005 


r 


i2q+ 020 

U.1Z»_ 020 


fl 1OO+0 020 
U.1ZO_ 020 


0.129t° ;gg 


n ioa+ 018 




25q5+ 0053 

U.zoyj_ 0047 


o 2 c ;qq +0 ' 0009 

U.zo»»_ 0009 


o.26oo±g:Sgg 






7qqo+ 0046 

U. I »£»->_0.0046 


7qq2+ 0026 

u. 1 S3i_ .0026 


7qq2+ 0027 
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on the last scattering surface, rescaled by a factor 100. 
The nuisance parameters are marginalized over in the fi- 
nal result. The marginalized 68.3% confidence level con- 
straints on cosmological parameters for Planck forecast 
only, Planck and LSS optimistic forecast, and Planck and 
LSS pessimistic forecast are listed in the second, third, 
and fourth columns of Tabic IIII1 respectively For com- 
parison we list the error bars obtained from a Fisher ma- 
trix analysis in the fifth column. 

The forecasted constraints on the plane n s -a s are 
shown on the top panel of Fig. [TJ Even in the pessimistic 
case, LSS surveys can improve the figurc-of-merit for 
{n s ,a s } (defined as FOM = 1/^/det Cov(n s , a s ) where 
Cov denotes the covariance matrix of the two parameters 
(63| ) by a factor of 2.2. Because the bias is unknown, 
large-scale structure data do not directly measure A s or 
us. However, our straw- man LSS survey can measure fi m 
to a much better accuracy, which can break the degen- 
eracy between f2 m and as that one typically finds using 
CMB data alone. This is shown in the bottom panel of 
Fig. [TJ We have checked that for Planck data r is almost 
orthogonal to n s and a s ; hence, our result is not sensitive 
to the fiducial value of r. 



C. Potential reconstruction 

Typically, slow-roll models predict a very small run- 
ning of the scalar spectral index, \a s \ < 10~ 3 . For these 
models the parametrization (jllj) is good enough. How- 
ever, this prediction is based on the assumption that the 
inflaton potential is "globally simple", i.e. that it does 
not abruptly change during the 60 e-folds of inflationary 
expansion. This assumption cannot always be tested, 
because the window of observed cosmological perturba- 
tions that ranges from the current Hubble radius to the 
scale where perturbations become non-linear is generated 
in about 10 e-folds of inflation. Perturbations on smaller 
scales, generated in the remaining ~ 50 e-folds of infla- 
tion, are difficult to study because of the complicated 
non-linear dynamics in the late-time universe. As a con- 



0.01 



-0.01 



Planck + LSS 
Planck 



68%, 95% CLs 




0.96 0.965 0.97 0.975 



0.81 



o s 0.8 



0.79 



Planck + LSS 
Planck 




, 95% CLs 



0.25 



0.26 



0.27 



FIG. 1. Marginalized forecasted posterior contours (68.3% 
and 95.4% CL) on the n 3 -a 3 plane (top panel) and Q,m-ag 
plane (bottom panel) obtained for Planck alone (blue-dashed 
lines) and combining Planck with a pessimistic LSS survey 
(red-filled contours). The white dots correspond to the fidu- 
cial model. 
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sequence, only a small piece of the potential is directly 
observable with the CMB and the large-scale structure. 

One way to constrain potentials independently of the 
slow-roll assumption is to compute V((f>) for a given 
model and constrain the parameters of this model. We 
will follow this approach in the next section. In this 
section we consider a complementary procedure, called 
potential reconstruction, which consists in studying only 
the part of the potential responsible for the observable 
window and relaxing the assumption of slow-roll outside 
this window. The remaining part of the potential can 
be of any shape, as long as it gives the right number of 
e-folds of expansion. 

Many methods have bee n prop osed to reconstruct the 
observable piece oiV(4>) jUaSil. All of them require 
some assumption regarding the smoothness of the poten- 
tial in the observable window, without which one could 
introduce arbitrarily many mass scales to build arbitrary 
shapes of the potential. The assumptions most used are: 
cutting at finite order some expansion in slow-roll param- 

i i, m, H3, m m 



ctcrs 



72H751 77, \Z 



between a finite number of bands 



S7R8 



interpolating 

~ Ml H3, m , 



direct inversion of the d spectrum [§1, [8l|-(86|] , or cutting 
off at finite order some expansion of the inflaton poten- 
tial V{4>) (6(1 [Z9|, [9H . Unfortunately, the results obtained 
depend on these assumptions and, for forecasts, also on 
the input fiducial parameters. 

Here we take the inflaton potential to be roughly 
smooth over scales of order ~ M p , i.e., 



Ml 



d n In V 



< 1 



n= 1,2,3. 



In the spirit of Refs. [HI, [T^] we assume that the first 
two slow-roll parameters are small. In particular ty = 
0(A) and r]v = 0(A) where A is a small parameter which 
observations - i.e. n s — 1 < 0.1 and r/16 < 0.1 - suggest 
is of order 0.1 or smaller. Note that we have not assumed 
that the running or higher-order slow-roll parameters are 
small. However, their smallness will be a consequence of 
our analysis. 

Before proceeding, let us first estimate the excursion 
of <p, which we call 5(f), during the 9 e-folds corresponding 
to the observable window 10 _4 Mpc _1 < aH < Mpc -1 . 
Using the slow-roll conditions, 5cj) is related to the varia- 
tion in the number of e-folds SN by [95| 



p V6.9 



1/2 



SN 



(22) 



Using SN — 9 we obtain 8<f> = 3A^/rM p , which implies 
that for reasonable values of r, the potential remains 
smooth inside all the observable window. 

We parameterize the observable piece of the potential 
by Taylor-expanding it in power series around the value 



of 



for which aH 



0.05 Mpc , and for 



simplicity we take 



v(<f>) 



E 

n=0 



0=0 



0. Thus, 



M r . 



M n 

P 



d n V 



(23) 

Without loss of generality, we can take V\ < and hence 
d<fi/dt > 0. With these definitions our smoothness prior 
translates into |U n /U| < O(l) for each n. 

Naively, one could truncate the expansion above at 
some arbitrary order n max and use Vq, V±, . . . , U„ max as 
parameters to constrain. However, it is important to un- 
derstand how these parameters enter into the expression 
of the power spectrum, the spectral index and running, 
which arc the quantities that are to be observed. Us- 
ing the smallness of ey and r?y , the inflationary power 
spectrum can be written as [18|, fl9j 



V 3 



12n 2 M^ 



391 - ~ 



(24) 



p=i 



V 



p-1 



o(\ 2 



where the potential V and its derivatives V n are com- 
puted at Hubble crossing. Here q p are known numerical 
coefficients of order unity [96|. Moreover, the spectral 
index is given by 



1-3 



M\ 9 V^ (Vx 



(21) while the running is 



V 



v, 



V, 



p+2 

V 



+ 0(X 2 ) , (25) 



p+3 

V 



+ 0(X 2 ) . 



(26) 



(Note that setting q p = for p > we recover the stan- 
dard slow- roll results of section IIII Al ) 

Given the assumption of smoothness introduced above 
and given that (Vi/V) 2 = 0(A), terms of order ~ 
T/(l/ ) 3 or higher are negligible in eqs. ((24)) - ([26)) and can 
be discarded. In other words, because in slow-roll infla- 
tion the excursion of 4>/M p during 10 e-folds is small, a 
local Taylor expansion truncated at 4th-order will be ac- 
curate enough to reproduce the correct power spectrum 
V s - The higher-order coefficients in the Taylor expansion 
are not measurable, unless they are much larger than 
unity, and they will not be included in our paramctriza- 
tion. 

Thus, in order to derive constraints on the potential 
we chose parameters that are closely related to those of 
the familiar CosmoMC, i.e. lnA s , r, n s and a s . These 
are 

In A, = In — ; — 77 , 

l2ir 2 M£V 2 ' 

2 t^2t 



n s = l + 2--3(- 



2qi 



V 2 



+ 2q; 



V{V4 
1/3 



(27) 
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where the right-hand sides of these expressions are eval- 
uated at 4> = 0, and 



'V V 



2q x 



V 



(28) 



0=0 



with 9 i = 1.063 and q 2 = 0.209 [jjj. A tilde is used to 
distinguish the parameters in eqs. (|2"7|) and (|2"8")l , defined 
in terms of the potential, from the observed quantities 
A s , n s and a s of the power spectrum in eq. and 
the tensor- to-scalar ratio r in eq. (|T3"|) . Finally, we still 
miss one independent parameter which we take to be 
proportional to V4/V and normalize in such a way that 
it is approximately equal to V4/V when r is about 0.1, 



= Va _f 
C4 = v 7 OJ 



(29) 



4=0 



We calculate V s and V t for each set of parameters 
in eqs. (j2"T)) - (j2"9"|) by numerically solving for the evolu- 
tion of scalar and tensor perturbations and hk for 
10 3 wavcnumbers log-uniformly distributed in —9 < 
ln[fc/Mpc _1 ] < 1 and by interpolating in between [§J 
Those parameters producing a non-monotonic potential 
or a non-inflationary phase (a < 0) are discarded. In 
particular, we study two examples with fiducial tensor- 
to-scalar ratio f = (top panel) and fiducial f = 0.128 
(bottom panel of Fig. [5]). The 68.3% CL constraints on 
the parameters in eqs. (|2"T]) - ([2T)| are reported in Table ITVl 
(In the fiducial r = case we report the 68.3% CL and 
95.4% CL upper bounds on f.) The uncertainties on 
these parameters - and on those not reported in this 
table - depend little on the fiducial value of r and are 
slightly worse than those obtained in section UlIBI in the 
case of slow-roll inflation (sec Table llllj) . which is ex- 
pected because we have introduced a new parameter, C4. 
Thus, the improvement upon the case of Planck alone 
is again mild. Moreover, the new parameter C4 will not 
be well constrained by future data. In conclusion, with 
Planck and LSS surveys we will be able to measure V3/V 

1 /2 

with an accuracy of ~ 0.1 (— ) and V4/V with an ac- 
curacy of ~ 0.5 (— J. Finally, for comparison with the 
slow-roll case, in Fig. [5] we show the constraints on the 
plane n s -a s . 

A final remark is in order here. Apart from the as- 
sumptions of smoothness discussed at the beginning of 
this subsection, here we also assumed uniform priors on 
In A, f, h s , a s and C4. For a well constrained parame- 
ter the choice of prior does not significantly change the 
results. This is the case, for instance, for the parameters 



hi As 



and a s used in our parametrization, which 



are closely related to observables. Note, however, that in 
other approaches to potential reconstruction, such as the 
flow-equation approach [64j , the cutoffs of the expansion 
and the priors/boundaries of parameters are often chosen 
arbitrarily. 



0.02 [ 
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Planck 

Planck + LSS 



-0.01 




0.96 
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0.02 [ 



0.01 . / 
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0.96 
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FIG. 2. 68.3% and 95% CL contours on the plane n s -a 3 
for Planck alone (blue dashed lines) and for Planck + LSS 
pessimistic forecast (filled red regions). On the top panel a 
fiducial f = is used, while in the bottom panel the fiducial 
value of f is 0.128 (such as for quadratic inflation). 



IV. CONSTRAINING FEATURES IN THE 
POTENTIAL 



In this section we study inflationary models that pre- 
dict features in the primordial scalar power spectrum. 
The amplitude of these features is small compared to the 
overall amplitude of the spectrum, typically of a few per- 
cent or less. Before dealing with concrete models, let us 
estimate the typical width in momentum space that can 
be measured by a perfect large-scale structure and CMB 
experiment. A similar discussion can be found in [36| . 
We will confirm these expectations more precisely in the 
following subsection when studying concrete models. 
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TABLE IV. Constraints on the inflationary parameters 



fiducial r = fiducial r = fiducial r = 0.128 fiducial r = 0.128 

Planck Planck + LSS pess. Planck Planck + LSS pess. 
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A. Analytical estimates 

1. Large scale structure survey 

For large-scale structure, the main limitation comes 
from the finite size of the survey. Indeed, the fact that 
the survey volume is finite sets a coherence wavelength 
fcmin ~ where V is the volume of the survey, 

below which all features get smeared [98[. Practically, 
the density field gets smoothed by a window function 
W determined by the survey geometry; in Fourier space 
this smoothing window is typically a compact ball with 
width of order fc m i n . For simplicity, let us neglect in this 
discussion the effect of shot noise and assume that the 
smoothing window is simply a Gaussian function with 
half width aw- In this case, for each wavemode k the 
estimated power spectrum is given by 



where 



P(k) 



\W(k)\' 



d 3 k' 



P{k')\W(\k-k'\ 



2tt(Tw 



exp 



(30) 



(31) 



of the bump is simply given by S In A s > N, 



-1/2 



where 



47rfc|£r B /fc3, 



is the number 



N B = Ank%a B V/(2Tr) 
of modes inside a shell of radius k B and width a B . By 
defining the relative width of the bump & \vtk B = a B /k B , 
this relation can be rewritten as 



Shi A s - ^/S In k B > 



kr\ 



3/2 



> 10' 



(33) 



where in the last inequality we have used the typical co- 
herence frequency of galaxy surveys such as the one con- 
sidered here, i.e. fc min = 10~ 3 /i/Mpc and k B < fc max = 
0.1/i/Mpc. In the opposite case, i.e. when the bump is 
much smaller than the coherence length, a B -C <Jw, one 
obtains 



P{k) ~ P (k) 



v B 



1 + 5 In A s ■ exp 



aw 



(k - k B f 



(34) 

Thus, the amplitude of the bump in the estimated power 
spectrum is suppressed by a factor <r B /aw while its 
width gets smeared out to aw- Although the bump 
cannot be resolved, it can be measured if ShiA s > 
(<Jw/<?b) 1/2 N b 1/2 » N B 1/2 . Using a B /k B = Slnk B , 
we can rewrite this relation as 



We take the half width of the Gaussian to be exactly 
equal to the coherence wavelength fc m ; n , aw = &min- 

Let us study the effect of the window function on the 
estimator in the case of a power spectrum with features. 
As a warm up, we first consider a Gaussian bump of 
relative amplitude ShiA s , centered at k = k B and of 
half width a B , i.e. 



P(k) = P (k) 



1 + S In A s ■ exp 



o 2 



2ai 



, (32) 



where Po(k) is a smooth power spectrum. The estimated 
power spectrum can be computed by plugging this ex- 
pression into eq. pop . In the limit where the smooth 
power spectrum varies little in an interval aw or a Bl 
one can easily compute the estimated power spectrum 
analytically and the calculation can be simplified by con- 
sidering wavemodes k 3> a B . As expected, for a window 
function much smaller than that of the bump, aw <C 
the estimator in eq. (|32[) reproduces the original power 
spectrum (|32[) . In this case the condition of observability 



Shi A. ■ 5 hik B > 



> 10" 



(35) 



Consider now the case of a power spectrum with peri- 
odic oscillations of frequency a B , i.e. 



P(k) = P (k) [1 + Sin A s ■ cos(fc/er B )] 



(36) 



Using eqs. pU]) and (f3Tj) we can compute the estimated 
power spectrum similarly to what done above. For fre- 
quency a B 3> aw one recovers eq. (|36p . while in the 
opposite case, a B <C aw, one obtains 



P(k)=P (k) 



1 + S In A s ■ exp 



-ITT ) cos(k/a B ) 

L<J B , 



(37) 

The frequency and phase of oscillations do not change 
but their amplitude gets exponentially suppressed for 
a B <C aw- The condition for observability is now ob- 
tained by comparing the estimated amplitude of the os- 



cillations with N, 



-1/2 



where AT to t is the total number of 
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modes in the survey, i.e. N to t ~ (k max /k m - m ) 3 . This then 
becomes 



5 In A s ■ exp 



2a% 



> 



3/2 



10" 



(38) 



In practice, to implement in our analysis the fact that 
there is a fundamental coherence wavemodc below which 
we cannot resolve, we have applied the following pre- 
scription. For each redshift bin we have sampled the 
wavenumber k in units of fc m i n (fc m ; n (z) = 2tt/V 1 ^ 3 (z)). 
For each (discrete) value of k we have computed the 
power spectrum using eqs. ([3"0")) and ([3T1) . where the width 
of the Gaussian window aw has been chosen to be 



V2hi2 /4tt\ 1/3 , 



0.302 K 



(39) 



In such a way, the real-space representation of the win- 
dow function, if cut off at its half-height, contains the 
same volume as that of the redshift bin. The fact that 
aw is smaller than k m - m allows us to neglect the overlap 
between window functions centered around neighboring 
values of k but slightly weakens the constraints discussed 
above. 

To summarize, using LSS data we will be able to re- 
solve features larger than roughly the coherence wave- 
mode of the survey, k mm , corresponding to a relative 
width <51nfc > 10 -2 , whose amplitude is at least of order 

(^min/^max) ^ 10 



2. Cosmic Microwave Background 

The coherence wavemode of the CMB, i.e. the smallest 
k mode probed by the temperature anisotropics, roughly 
corresponds to k mm ~ D~J., where D I0C is the comoving 
distance to the last-scattering surface. More specifically, 
D rcc = To — T rC c, where To and r rec are the conformal time 
today and at recombination, respectively. Thus, we ex- 
pect that features in the primordial power spectrum with 
a width in momentum space smaller than D~J, cannot be 
measured. This implies that we can only resolve features 
with relative width such that Sink > k m i n /k max ~ f^, 
where £ max is the largest multipole of the CMB map. 
However, this naive estimate is too optimistic. Indeed, 
the projection from momentum space to multipole space 
degrades the amplitude of features, hampering their mea- 
surement. The order of magnitude of this effect can be 
schematically computed by the following example. 

To simplify things as much as possible, let us consider 
the angular spectrum of the temperature anisotropies in 
the Sachs- Wolfe, instantaneous recombination and flat- 
sky limit. This is given by 



1 



187T£>2 



-foo 



dyP(k=^{ 2 /D? cc +y*). (4Q) 



Plugging the expression for the power spectrum with os- 
cillations, cq. (|36[) . into this equation, one can solve the 



integral on the right-hand side for a scale invariant spec- 
tral index of primordial fluctuations. Its solution can be 
written in terms of a Mcijcr G-function 199], but it is 
more enlightening to use its asymptotic expansion in the 
large x = £/aBD ICC limit, by which we can recast eq. (|40|) 
as 38] 



1/2 



1 + 5\nA : 



x(t) + 



(41) 

where Co.e is the angular CMB spectrum for a smooth 
initial power spectrum. Thus, apart from a phase shift, 
the oscillatory behavior does not change. However, the 
amplitude of the oscillations gets suppressed by a fac- 
tor ~ ((Tg-Drec/f) 1 / 2 [38|. For oscillations with constant 
frequency in k the suppression factor is proportional to 
£~ x / 2 and hence becomes more important at large £. Note 
however that for oscillations with constant frequency in 
lnfc, such as in axion monodromy, the suppression factor 
is independent of £ and, using as — (£/D lec )Shik, it is 
simply given by \J 6 In k. 

To measure these oscillations, their relative amplitude 
in the Ci must be larger than the inverse square-root of 
the number of pixels. From the discussion above this is 
given by 



SlnAs •V<5hifc>5-1(T 



(42) 



where we have used Z max = 1500 (for larger £'s, the noise 
is larger than the signal). 

As shown in Ref. [38(, there is also a smearing effect 
due to gravitational lensing starting at £ ~ 10 3 . As most 
of our constraints come from £ < 10 3 , where CMB lensing 
is negligible, we did not include this effect in our analysis. 

For the analysis that will be discussed in the rest of 
this section, the presence of small features and oscilla- 
tions makes the calculation of the angular power spec- 
trum challenging. Indeed, in order to speed up the com- 
putation, in publicly available CMB codes such as CAMB 
IH, CLASS mEl] or CMBEasy [H], the Q's arc com- 
puted by sampling the primordial power spectrum with 
steps longer than the width of these features and using a 
certain number of approximations which in our case are 
not accurate enough. For this reason, one of the authors 
of this article (ZH) has developed a new CMB Boltzmann 
code which is adapted to treat very sharp features in the 
primordial power spectrum at a reasonable computing 
time. This code is written in the Newtonian gauge and 
agrees well with other codes such as CAMB and CLASS. 
The code is presented in [46| but we outline its main 
features in the appendix, Sec. 



B. Ringing features 

As an example of model of inflation predicting a sharp 
feature in the power spectrum we consider the model 
first presented by Starobinsky in [2(J. The feature is 
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produced by a discontinuity in the first derivative of a 
linear potential. More precisely, this can be defined by 
its amplitude at <f> — 0* (the value of the field at Hubblc- 
exit of the pivot scale fc*), i.e. V* = V (</>*), and its first 
derivative given by 



M, 



dV = \ Vi,. 



if, 
if, 



> </>± 



(43) 



where Vq : _ and Vi + arc constants. A closed form of 
the potential can be obtained by integrating dV/dcj) from 
= 0*. With this definition the inflaton mass d 2 V/d<fi 2 
diverges at <fi = <f>±, which is theoretically fine but 
numerically problematic. Thus, in practice we take a 
slightly modified potential which smoothly interpolates 
the jump of the first derivative: we use eq. (|43|) for 
- > eM p and we linearly interpolate 

— < eM v , where e is a small 



dV/d(j) for \<t>- 4>± 



Vi,-+Vi, 



parameter. In the example below we consider e = 0.01. 
We have verified that V s and Vt are stable to changes of 
s around this value and that they do not depend on the 
choice of e as long as e < 1. 

We constrain the following parameters, constructed 
from V*, Vi -, Via and cj)±, 



In A, 



In 



V? 



6n 2 M^(V^ 
Vi,- 



V, 



vl + ) 



V. 



which approximate the average of A s and n s 
tively, and 



Sn s 



Vi,- 
V* 



Vi, 



V, 



ln(fc ring /fc») = 4>i 



1 - n s 



(44) 

(45) 
respec- 

(46) 
(47) 



which parameterize the variation of the spectral index 
and the position where the ringing feature starts. As 
above, the tilde is used to distinguish quantities defined 
by the potential from those defined directly in terms 
of observables. For each set of parameters given in 
eqs. (|3i )) -(|3? |) we solve V*, Vi-, Vi i+ and <f>± and cal- 
culate scalar and tensor perturbations and hk for 
10 3 wavenumbers log-uniformly distributed in —9 < 
ln^/Mpc" 1 ] < 1. 

As an example, in Fig. [3] we show the scalar power 
spectrum calculated by numerically solving the evolu- 
tion equations for scalar perturbations. The solid red 
line corresponds to the case \nA s = 3.02, n s = 0.975, 
Sn a = 0.002 and ln(fc r i ng /fc*) = —2, which will be stud- 
ied below. The slow-roll approximation [59| is also shown 
for comparison and as a check of the numerical accuracy. 
To see the effect of varying Sn s on the power spectrum, 
in the same plot we also show the case for Sn s = 0.004 in 
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- numeric, 8 n s = 0.002 
slowroll, 5 n s = 0.002 
numeric, 8 n = 0.004 
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FIG. 3. The primordial power spectrum lnVs, computed by 
numerically solving the scalar perturbation equations, for a 
step in dV/d(j> in the Starobinsky model. The solid red line 
is the power spectrum for the model with \nA s = 3.02, n a — 
0.975, Sn s = 0.002 and ln(fc ring /fc*) = -2; the dotted line is 
the slow-roll approximation for the same model. The dashed 
cyan line is the power spectrum computed for a different value 
of <5n s , i.e. Sn a = 0.004. 



cyan. Typically, in this model the relative width of the 
feature is of order 1. Its exact value is related to 6n s , 
but as long as we consider only small values of 5n s it is 
not too sensitive to it. Moreover, using eqs. (|4"4")) - (|4l)|) . 
we can relate the relative jump in the amplitude of the 
power spectrum to 5n s , i.e. 



5 In A, = 



SAs 
A s 



25ns 
1 — h 



In the case of the example above one finds S In A s 
which is confirmed by the numerical calculation. 



(48) 



0.16, 



1. Null test 

Let us first find the upper bound on 5n a that can be ob- 
tained in absence of a ringing feature, i.e. by choosing the 
fiducial 5n s = 0. For this null test we produce mock data 
using a fiducial linear potential model with ln(10 10 A s ) = 
3.02 and n s = 0.975. As a prior on m(fc r i n g/fc*) we take 
—4 < ln(fc r i ng /fc*) < 2, which roughly corresponds to the 
range of scales probed by Planck and the LSS galaxy sur- 
vey, i.e. 10" 3 Mpc -1 < k < 0.37 Mpc -1 . Note that in 
the null-test case one has to chose this prior with care. 
Indeed, with a bad choice of priors the posterior proba- 
bility may be dominated by a fraction of the volume of 
the parameter space which is beyond the observational 
window and thus insensitive to the data. 

We show the constraints using CMB + LSS in Table IVl 
and the ln(fc r i ng /fc Ht )-5n s plane in Fig. 01 The standard 
cosmological parameters such as Qbh 2 , Q c h 2 , 9 and r re 
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TABLE V. Null-test on Starobinsky model 
Parameter Planck + LSS pess. 



Q b h 2 

n c h 2 



0920+ 0001 
u.uzzu_ 0001 

n i 1 280+ ' 000 ' 



1.0462 



+0.0002 



o.o9otg-°» 



ln(fc r i ng /fc*) unconstrained 

lnClO 10 ^) 3.013ig"gS 

Sn s ^ 0000 +o.ooo25+o.ooo66 



0.001 



0.0005 




Planck + LSS 



-4 -3-2-10 1 2 
In (k ring /k t ) 

FIG. 4. The marginalized 68.3% CL and 95.4% CL con- 
straints on Sn s as a function of ln(fc r mg/fc*) for ringing features 
in the Starobinsky model. The mock data are produced using 
a fiducial linear potential model with ln(10 10 J 4 s ), n s = 0.975, 
and no feature (Sn 3 — 0). The prior on ln(fc r i ng /fc*) is 
-4 < ln(fc ri n g /fc*) < 2. 



are not affected by the ringing feature, so that choosing a 
model with an incorrect value of 8n s or ln(fc r jng/&*) has 
almost no impact on these parameters. The upper bound 
of 6n s depends on m(fc r i ng /A;*). Indeed, since on smaller 
scales more Fourier modes are probed by CMB and large- 
scale structure, the constraint on Sn s is better for larger 
ln(fc r i ng /fc*), as shown in Fig. |4] It is necessary to use 
both CMB and large-scale structure to cover the wide 
range of scales that we have chosen in the prior. Using 
only one of them would lead to posteriors dominated by 
an unconstrained region in the parameter space. Thus, in 
this case we do not compute the CMB-only constraints. 



2. Measuring the ringing feature 

Now we study how well CMB and large-scale structure 
can detect a ringing feature and as fiducial values we 
choose 8n s — 0.002 and ln(fc r jng/&*) = —2, corresponding 



TABLE VI. Constraints on Starobinsky model 
Planck only Planck + LSS pess. 
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ln(fcring/fc») 

fn(10 10 i s ) 
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-2 nnn +0 01 ° 
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o n21 +0.011 
O.U^±_q 01() 
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.o nn4+ - 007 



3.022 



+ 0.00K 



0.00200±™ 0.00203±°= 



0.003 



0.002 



0.001 




-2.02 



In (k rme /k«) 



FIG. 5. The 68.3% CL and 95.4% CL constraints on Sn s 
and ln(fc r i ng /fc») with Planck-only (blue) and Planck plus LSS 
survey pessimistic forecast (red), for ringing features in the 
Starobinsky model. The white point is the fiducial model. 



to the case plotted in Fig. [3] The fiducial values on the 
other parameters are chosen as in the null-test case. 

The forecasted constraints on the cosmological param- 
eters are listed in Table IVI1 The constraints on the pa- 
rameters not listed in this table are very close to those 
for the null case, listed in Table [V] The constraints on 
the ln(/c r i ng /fc*)-(5n s plane are shown in Fig. [SJ As we 
explained in the previous subsection, when mapping the 
power spectrum to the angular Ce, the narrowest features 
in the power spectrum get suppressed. Thus, the CMB 
experiment is likely to capture only the first two or three 
oscillations of the ringing features while a galaxy survey 
captures more oscillations. Indeed, the LSS data improve 
the total FOM for 8n s , In(fc r i ng /A;*), n s and ln^4 s by a 
factor of about 5. Note that the constraints for LSS can 
be qualitatively derived using eq. ([53")) . Indeed, using that 
the relative amplitude of the features is of order 1, this 
equation and eq. (|48[) imply that the typical constraint 
on 8?i s is of order a(Sn s ) ~ (& m in/^ring) 3 / 2 (l — n s ). By 
plugging the numbers in we find a(Sn s ) ~ 0.0015, which 
confirms the numerical results. 
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C. Constraining oscillations: axion monodromy 

For the axion monodromy model we use the parame- 
terized potential given in Ref. (40j . 



v{<j>) = (j? 



4> + bf cos 



(49) 



where p, parameterizes the scale of inflation, 6 is a small 
parameter (monotonicity of the potential requires b <C 1), 
/ is a - typically sub-Planckian - frequency, / <C M p , and 
inflation requires <f> 3> M p . In this model inflation ends 
when <fi ~ 0. As usual, <j> t denotes the value of the field 
when the pivot scale exits the Hubble radius, which for 



60 e-folds of inflation is 



11M„ 



Let us now map a variation of / in field space, 5(j) = 
/, into the corresponding variation in momentum space. 
We work at lowest order in the slow-roll approximation, 
where k ~ e m H with H constant. In this case one finds 



5\nk = -^c 



if- 



(50) 



Furthermore, neglecting the oscillations of the potential 
in the Klein-Gordon equation of the field <fi one obtains 
3Hip = —p 3 and replacing H using the Friedmann equa- 
tion (3M%H 2 = p 3 (j>) one gets H/<j> = -cp/M 2 . With 0* 
as typical value for the scalar field <j> we obtain the width 
of oscillation of the power spectrum in In k, 



Sin k 



f<t>* 
M 2 



(51) 



which will be used as a parameter to constrain. 

We parameterize the primordial power spectrum using 
the analytical approximation given in Ref. |40j , 



V s {k) = A s ( A 



1 + Sn s cos ( — + Atp 



(52) 



where 



TABLE VII. Null test on axion monodromy 



Parameter Planck only 



Planck + LSS pess. 



Q c h 2 



ln(10 10 A 3 



2 2 00+ 00011 

u - uzzuu -0. 00011 



0.021996 4 



0.0009 
+0.00021 
0.00022 
+0.004 
0.004 



1.04618 
0.090 

q no7 +0.008 
°- uz ' -0.008 



1.046181™ 

0-0901°:™ 

3.0268±g;ggS? 



9749+ 0023 

00 +0 0046+0 ' 0132 



0.9751 
0.0000 



+0.0018 
-0.0018 

+0.0015+0.0046 



becomes 



V s {k) = A s [ — 



1 + Sn s cos 



lnfc/fc* 



Sink 

' (55) 

where the phase tp is given by tp = — Atp— <f>*/ f. We have 
verified that Eq. ((55|) . which does not assume the slow- 
roll approximation, for b -C 1 and dink > 0.003 agrees 
to the percent level with the power spectrum obtained 
by numerically solving the scalar perturbation equations 
during inflation. 

In the following we will constrain 5 parameters. Two 
of them, the width of the oscillations Sink, defined in 
cq. ([5TI)) . and the free phase <p, have been already intro- 
duced. We allow the phase if to vary between — it and 
7r. We will then constrain the amplitude and spectral in- 
dex of the oscillations, which are defined in terms of the 
physical parameters </>* and /i 3 as 



,3^3 



12tt 2 M6 



= 1 



3M p 2 



(56) 



and the deviation from the spectral index n s , which is 
related to the physical parameter b and / in the poten- 
tial by cq. (f53|) . For the MCMC calculation we assume 
uniform priors on these five parameters. For the tensor 
spectrum we use that, at lowest-order in the slow-roll 
approximation, r = (8/3)(l — n s ) and n t = — r/8. 



Sn, 



126 



1 



'^cothf^ 
8 V2/0* 



f<f>* 
M 2 



(53) 



The phase Aip originates from the uncertainty both in 
the exact number of e-fold from the end of inflation and, 
on a more microscopical level, in the exact sinusoidal 
modulation of the scalar potential arising in the string 
theory construction. 

By integrating the first equality of eq. ([50| using again 
H/tj> = -(f>/M 2 one obtains [H 



4>k = V <j>% - 2 In k/K — 4 1 * 



In fc/fc* 



(54) 



which can be plugged into the argument of the cosine in 
cq. (j52"j) . Rewriting this argument using cq. ([5"Tj) . cq. (f5"2")) 



1. Null test 

As before, we start by considering the null test. Mock 
data is generated using a fiducial power-law spectrum 
with n s = 0.975 and ln(10 10 i s ) = 3.027. The constraints 
on the parameters are given in Table IVIII As the fre- 
quency Sink and phase tp remain unconstrained, they 
are not shown in this table. The constraints in the plane 
S In k-Sn s are shown in Fig. [6] As expected, in the Planck- 
only case the limit on Sn s depends on the frequency Sink. 
For Sink = 0.1 we find Sn s < 0.2 at 68.3% CL, which 
is about what expected from eq. (|4"2"j) . This constraint 
weakens at smaller S In k. Adding a large-scale structure 
survey with specifications such as those considered here 
improves these constraints, which for S In k > 0.01 are 
independent of Sink, as expected from the discussion of 
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o 
to 



0.02 



0.01 




FIG. 6. The marginalized posterior 68.3%CL and 95.4% 
CL constraints on Sn„ and S In k for Planck-only (blue) and 
Planck plus a LSS survey (red). Note that there are some 
unavoidable resolution effects in this representation of the 
constraints: the 1 and 2-cr contour lines should go infinity 
for tHnfc — > 0, as for very small values of <51nfc one does not 
have any constraints on oscillations. The fact that contour 
lines remain finite is due to a practical choice of finite range 
for <51nfc: our prior on Sink is 0.001 < Sink < 0.1 (for 8n a 
we have taken < Sn 3 < 0.05). This problem disappears, of 
course, in the case of a detection. 



0.02 



Sec. IIV A 11 Indeed, when Sink becomes smaller than 
fcmin/^max, which roughly corresponds to the relative co- 
herence length of the survey, the amplitude of oscillations 
of the estimated power spectrum remains unconstrained. 

Note that one should be aware of prior-dependence 
when interpreting the contours shown in Fig. [6] For in- 
stance, one would expect to see the 1 and 2-cr contour 
lines to go to infinity for 8 In k — > 0, as for very small 
values of 5 In k one does not have any constraints on os- 
cillations. The fact that contour lines remain finite is due 
to a practical choice of finite range for Sink, i.e. 0.001 < 
8 Ink < 0.1 (for Sn s we have taken < 8n s < 0.05). This 
ambiguity disappears, of course, in the case of detection 
(see below). 



2. Measuring the oscillations 

Let us now study how well we can constrain the ax- 
ion monodromy parameters in the presence of oscillations 
with amplitude 8n s = 0.01 and phase ip — 0. The choice 
of the fiducial value of ip is irrelevant for this analysis. For 
the fiducial frequency we look at two cases: <51nfc = 0.1 
and 8 Ink = 0.01. The other fiducial parameters are as 
for the null-test case above. The constraints on all pa- 
rameters are collected in Table lynTl 

For 8 Ink = 0.1, Planck-only mock data give a de- 
tection of oscillations, which is what expected from the 
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FIG. 7. The 68.3% CL and 95.4% CL constraints on 5 Ink and 
8n s with CMB only (blue) or CMB+LSS forecast (red). The 
white points are the fiducial models: on the top panel this is 
Sink = 0.1 and 8n s = 0.01; on the bottom panel 5 lnfc = 0.01 
and 5n s = 0.01. 



naive estimate, eq. (|42[) . The inclusion of LSS data can 
significantly improve the constraints on <5 In k and 8n Sl as 
shown in Fig. [7] Note that, in the case of a detection, the 
inclusion and marginalization over the phase parameter 
ip is more important than in the null-test case. Indeed, as 
shown in Fig.[5]there is a strong correlation between 8 In k 
and ip, which is absent in the null-test case. The usual 
way of updating the proposal matrix using the covariance 
matrix, e.g., the method used in CosmoMC, fails due to 
the multiple branches of posterior shown in Fig. [5J To 
solve the problem, we have written a new MCMC engine 
that allows parameters to have periodic boundary condi- 
tions. The MCMC engine is combined with the new Ce 
integrator described in the appendix j4fj| . 

Let us turn to another choice of fiducial value for the 
frequency, Sink = 0.01. As shown in the bottom panel 
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TABLE VIII. Constraints on axion monodromy 
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FIG. 8. The forecasted 68.3% CL and 95.4% CL constraints 
on the (ip-8 In k) plane for Planck-only (blue) and Planck plus 
a LSS survey (red). We considered fiducial values of if = 
and 8 In k = 0.1. 



of Fig. in this case Planck data cannot measure the 
oscillations. This is due to the suppression of features in 
the Cg spectrum occurring in the projection from k to £ 
space. To explicitly sec this suppression, we compute the 
CMB temperature-temperature angular power spectrum 
In Cj T for a primordial power spectrum InPs with and 
without monodromy sinusoidal modulation. The differ- 
ences between the two, for Sink = 0.1, 0.01 and 0.001, 
are shown in Fig. [5] Apart from the modulation due 
to the acoustic transfer, which cannot be capture by the 
simple analytic argument of Sec. lIV A~2l(38l ]. the suppres- 
sion due to the projection from k to £ space shown on 
the right-hand panels is roughly given by s/S In k, which 
confirms our discussion of Sec. IIV A 21 Moreover, at low 
^ < (Sink) -1 , the oscillations in /c-space disappear 
in £-space: indeed, because of the discreteness of t they 
simply cannot be sampled. 

Note that the error bar on the parameter S In k is ex- 
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FIG. 9. The differences in lnV a (left panels) or lnCj T 
(right panels) between a fiducial axion monodromy model 
with In (l0 10 yl s ) = 3.027, n s = 0.975, amplitude of cosine 
modulation Sn s = 0.01, phase <p — and a smooth power-law 
spectrum with the same A s and n a . From top to bottom, 
<51nfc = 0.1, 0.01, 0.001, respectively. 



tremely tiny in the case of detection from CMB + LSS 
shown in the lower panel of Fig. [JJ One can estimate this 
error as 



c<5 in fe ~ N os l 5h\k. 



(57) 
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FIG. 10. Slices of the CMB likelihood (dashed blue lines) and 
the large-scale structure likelihood (solid red lines) for varying 
Sink. The other parameters are all fixed at the fiducial val- 
ues. The fiducial Sink = 0.01 is where the spike of likelihood 
presents. Both likelihoods are normalized to 1 at their max- 
ima. The main frame shows the large-likelihood part that 
includes the CMB likelihood slice and a sharp spike of the 
large-scale structure one, while the lower-right panel zooms 
out the global structure of the large-scale structure likelihood 
slice. 



where N osc is the number of oscillations in the observable 
window in lnfc. This estimate can be found by consider- 
ing the shift in S In fc inducing a change in iV osc of order 
unity. Indeed, a much larger change of N osc makes most 
of the oscillations out of phase, hindering the detection. 
Thus, using that N osc oc (<51nfc) _1 we find eq. (|57|) . Both 
panels of Fig. [7] confirm this estimate. 

Note also that from eq. (j42|) we would expect a de- 
tection. However, this estimate is too naive to be valid 
also in the regime where this equation is only marginally 
satisfied, and one has to rely on the numerical analysis. 
Adding the LSS improves considerably the figure. In- 
deed, in this case we can detect oscillations and obtain 
rather tight constraints on both Sink and 8n s . 

If we vary 5 In fc and keep all the other parameters fixed 
at their fiducial values, the form of the likelihood as a 
function of 8 In fc is very different from a Gaussian: it is 
typically sharply spiked around the fiducial value of 8 In k 
and remains almost constant far from the fiducial value 
of 8 lnfc. An example is shown in Fig. QUI This form of 
likelihood cannot be studied by standard Fisher matrix 
analysis. Indeed, in this case searching the maximum of 
the likelihood is more similar to the process of tuning a 
radio channel to search for the correct frequency. Since a 
tiny change of 8 In fc can destroy the tuned-in resonance 
between the data and the theoretical prediction, the error 
bar on <51nfc is either tiny in the case of a detection, or 
very large in the case of no detection. 

This is exactly what happens in the case of fiducial 
8 In fc = 0.01 discussed above. In the Planck-only case the 



parameter space is dominated by the tuned-out situation, 
where the frequency and phase of the oscillations due to 
monodromy in Cf Y and C* Y of eq. ([T]) arc different. 
In this case the % 2 of (HJ is minimized by models with 
small 8n s (\ 2 will be roughly proportional to (8n s ) 2 + 
(8n s ) 2 , where the cross term vanishes due to uncorrelatcd 
phases). Thus, after marginalization over all the other 
parameters the fiducial value Sn s = 0.01 is not a better 
fit than 8n s = 0. Indeed, as shown in Fig. [7j the fiducial 
model is not even in the marginalized 1-er contour. 



V. CONCLUSION 

In this paper we studied how planned stage IV galaxy 
rcdshift surveys (according to the classification of fl5| ) 
with characteristics not too dissimilar from those of, e.g., 
the Euclid mission, will improve the constraints on infla- 
tionary models that one obtains from CMB data alone, in 
various situations. For slow-roll models and in the more 
general case of a smooth inflaton potential over the field 
range of a Planck mass (see Sec. . we have studied the 
forecasted constrains on the spectral index and its run- 
ning. From the expected signal-to-noise of the Planck 
satellite data we have found that by adding LSS galaxy 
survey data we obtained an improve ment of about a fac- 
tor of two, in broad agreement with |l00j . 

Then, in Sec. II VI we have considered models predicting 
sharp features in the power spectrum. In this case galaxy 
surveys will be crucial to detect and measure features. 
Indeed, what we measure with the CMB is the angular 
power spectrum of the anisotropics in the 2-D multipolc 
space, which is a projection of the power spectrum in 
the 3-D momentum space. As explained in Sec. IIVI fea- 
tures at large £'s and for small width in momentum space 
get smoothed during this projection (65J. Thus, future 
galaxy surveys will be able to measure features that are 
invisible to CMB observations. The main limitation on 
the width of features measured using LSS comes from 
the size of the volume of the survey, the smallest de- 
tectable feature being of the order of the inverse cubic 
root of this volume. The other limitation comes from the 
number of modes contained in this volume. In order to 
improve measurements of oscillations, it will be impor- 
tant to have a better control on the maximum wavemode 
that can be used in the analysis. This can be done only 
by reaching a better understanding on rcdshift distor- 
tions, non-linear physics at small scales and galaxy bias. 
These findings arc robust to small changes of the survey 
characteristic especially for models which give a smooth 
power spectrum. Models with sharp features or oscil- 
lations however are very sensitive to the fraction of the 
observed sky (i.e. to the survey volume) as it is clear from 
Sec. lIVAll 

For models predicting oscillations in the power spec- 
trum we have found that future LSS will push the ampli- 
tude of observability down to 8n s ~ 10 -3 , provided that 
the frequency of oscillations is 0.01 < 8 lnfc < 1. (For 
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5 In k ^> 1 the oscillations again become unmeasurable 
due to the limited range of observable scales.) Quantita- 
tively similar conclusions have also been reached earlier 
in [36[. 

Note that hints of monodromy modulation in WMAP 
and ACT CM B power spectrum has been recently re- 
ported in Refs. fnnlfio^ . Since the best-fit S In k claimed 



in Ref. 101 



is of order 10~ 2 , caution needs to be taken 
for the accuracy of Cg integration, the marginalization 
and the MCMC convergency. Most notably, a signifi- 
cant improvement of \ 2 a t some best-fit point does not 
necessarily indicate a detection of oscillations, as the like- 
lihood cannot be approximated as a multivariate Gaus- 
sian. Indeed, the absolute height of the likelihood spike 
alone cannot determine whether there is a detection or 
not. What matters is the product of the height of the 
likelihood spike and the volume that it occupies - see 
Fig. [TU] and the discussion in the last few paragraphs in 
subsection IIVC 21 

For this analysis one of us (ZH) has developed a Ci 
integrator, a modified MCMC engine that allows param- 
eters to have periodic boundary conditions and a forecast 
mock data generator that does the calculation described 
in Section [TTJ We briefly review their main features in 
the appendix. The package with all these codes is a self- 
contained tool to forecast constraints using large-scale 
structure, CMB and supernovae - although the latter 
have not been used in this paper. It is publicly released at 
|http : / / www . cita . ut oronto . c a/ ~zqhuang/ CosmoLib 
and described in more details in an accompanying paper 

El. 
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Appendix A: The Numerical Package 

The modified MCMC engine and the new integrator to 
compute more accurately the CMB angular power spec- 
tra for models with rapid oscillations is presented in de- 
tail in Ref. [46|. Here we summarize its main features. 

For the MCMC engine, to achieve quicker convergence 
in the case that the likelihood is a periodic function of 
some parameter ip, in each random-walk step we have 
mapped the proposed new value of p into [— tt, 7r) by 
adding integer number of 27r's to tp, instead of using the 
usual rule that tp is rejected when it exceeds the bound- 
aries. 

The challenging task is to accurately compute the an- 
gular CMB spectrum for each multipole £, in the presence 
of oscillations in the primordial power spectrum V s {k). 
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FIG. 11. The temperature transfer function A* for a fixed 
£ = 300 



The angular spectrum for the temperature anisotropics 
is given by 



A k e \ 2 V s (k)dlnk , 



(Al) 



where A| is the temperature transfer function, which is 
typically oscillatory. Since A^ is evaluated numerically, 
its sampling is time consuming. Indeed, in modern fast 
CMB codes - such as CAMB [II, CLASS 0, E|, or 
CMBEasy [45| - the integral above is computed by sam- 
pling A{? using a step size which can be typically much 
larger than the oscillation period. For instance, in Fig.QT] 
we show an example of A^ for a fixed £ = 300. A typical 
sampling scheme is shown by the red solid triangles in 
the upper-right panel, which zooms-in part of the figure. 
If V s (\~a.k) is a smooth function - i.e. a low-pass window 
- sparse sampling of A^ is enough. 

However, when V s {k) has local sharp features, at least 
a few samples of Aj? are needed per dink, the typical 
width of the features. For instance, if our goal is to 
sample features with width i51nfc ~ 10 -3 , the required 
sampling frequency is typically ~ 20 to 100 times higher 
than that used for a smooth V s (k). Furthermore, as we 
wish to compute the C/s for each £ rather than interpo- 
lating it over few tens of £'s, its computation is typically 
few thousands times more expensive than in the smooth- 
V s (k) case. A final complication is due to the fact that, 
if all the transfer functions and the precomputed je(x) 
tables are to be stored, one has also to face a memory 
barrier that cannot be easily bypassed. For these reasons, 
simply increasing the I and k resolution in standard codes 
such as CAMB, CLASS or CMBEasy, will not be efficient 
enough for the purpose of scanning the whole parameter 
space. 

Here we compute the C/s in equation (|Al[) £ by £. We 
chose the sampling period in In k to be at least a few times 
smaller than <51nfc. To speed up the computation and 
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avoid memory problems, our strategy is the following. 
We first compute two neighboring C/a by brute force 
(see below). Two arrays of spherical Bcssel functions 
jt+i [k(ro — t)] and jg [fc(ro — r)] are stored in the memory 
for each (/c,r) sample used in the computation. Then 
we compute Ct-i. To do that, we directly compute 
using the recurrence relation 

21 + 1 

3i-i(x) = ji{x) - jt + i{x) . (A2) 

x 

The values of jg+\ are then discarded to free the memory. 
Again, using ji and jt-\ we then calculate ji-2 and hence 
C(-2- This downward iteration is very stable for a few 



tens of steps, after which we need to recompute another 
couple of neighboring CV's and iterate downward again. 

The initial neighboring j/s are calculated using pre- 
computed 25-th order Chebyshev fitting formulas. (For 
the rapidly oscillating part at x 3> I, the phase and am- 
plitude of oscillations are fitted using Chebyshev poly- 
nomials.) Chebyshev fitting is slightly slower than the 
cubic-spline fitting used in other public CMB codes, but 
it is more memory-efficient and more accurate - it has an 
accuracy ~ 10 -8 - and allows more downward iterative 
steps. Finally, note that the algorithm proposed here is 
more efficient both CPU-wise and memory-wise, enhanc- 
ing the speed of of £-by-£ computation of Cgs by a factor 
of - 10 to 30. 
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